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Simulations of bent-double radio sources in galaxy groups 
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ABSTRACT 

Bent-double radio sources have been used as a probe to measure the density 
of intergalactic gas in galaxy groups. We carry out a series of high-resolution, 3- 
dimensional simulations of AGN jets moving through an external medium with a 
constant density in order to develop a general formula for the radius of curvature 
of the jets, and to determine how accurately the density of the intra-group medium 
(IGM) can be measured. Our simulations produce curved jets ending in bright radio 
lobes with an extended trail of low surface brightness radio emission. The radius of 
curvature of the jets varies with time by only about 25%. The radio trail seen in our 
simulations is typically not detected in known sources, but may be detectable in lower 
resolution radio observations. The length of this tail can be used to determine the 
age of the AGN. We also use our simulation data to derive a formula for the kinetic 
luminosity of observed jets in terms of the radius of curvature and jet pressure. In 
characterizing how well observations can measure the IGM density, we find that the 
limited resolution of typical radio observations leads to a systematic under-estimate 
of the IGM density of about 50%. The unknown angles between the observer and the 
direction of jet propagation and direction of AGN motion through the IGM leads to 
an uncertainty of about ±50% in estimates of the IGM density. Previous conclusions 
drawn using these sources, indicating that galaxy groups contain significant reservoirs 
of baryons in their IGM, are still valid when considering this level of uncertainty. 
In addition, we model the X-ray emission expected from bent-double radio sources. 
We find that known sources in reasonably dense environments should be detectable 
in ~ 100 ks Chandra observations. X-ray observations of these sources would place 
constraints on the IGM density and AGN velocity that are complementary to radio 
observations. 



1 INTRODUCTION 



Observations of the local universe support the conclusions 
drawn by the hierarchical theory of structure formation 
that the majority of galaxies reside in groups, dynam ically- 
boun d systems which span a wide range of properties l|Tullvl 
1 19871 ). The intra-group medium (IGM) contained within 
these groups likely contains a significant fraction of the bary- 
onic content of the universe |Fukugita et alj|l998h . Obser- 
vations measuring the baryon content of the local universe 
account for approximately a third of the baryon density ob- 
served at high redshift (z = 2 — 4 ) in the form of stars, cold 
gas, and hot X-ray -emitting gas (|Fukugita fc Peebles! 1200^ ; 
IStocke et al.l 120041 ). Numerical simulations predict a signif- 



icant fraction of the "missing" baryons, around 40% of the 
total baryonic content of the local universe, may be con- 
tained in the warm-hot intergalactic mediu m (WHIM) with 
a temperature range of T ~ 10 5 - 10 7 K (|Cen fc Ostrikerl 
l2006l ; lDave et al.ll200lf ). 

The primary observational tool for studying the WHIM 
is curre ntly absorption-line s pectroscopy of low-redshift 
quasars IjNaravanan et al J20ld ). These observations are lim- 
ited to groups falling along the line of sight of a sufficiently 
bright quasar, and density measurements depend on esti- 
mates of the systems' spatial extent, metallicity and ioniz- 
ing fraction, yielding total IGM de nsity measurement s on 
10" 4 - 10" 5 cm -3 (|Pisano et alj l2004h . X- 
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ray measurements are inherently limited to higher temper- 
ature groups, generally containing at least one early-type 
galaxy. A dynamical mass can be estimated by making an as- 
sumption about the geometrical distribution of X-ray emit- 
ting gas in the group. However, measurements of X-ray sur- 
face brightness drop below measurable levels well within the 
virial radius, requiring additional assumptions i n order to 
extra polate the total mass estimate of the group |Mulchaevl 
l2000h . 

Observations of bent-double radio sources in galaxy 
groups serve as a valuable method of measuring IGM den- 
sities, and recently this method has been used to find 
total IGM densities in groups of n ~ 2 x 10~ — 3 x 
10 -3 cm" 3 (|Freeland et alj|2008l . I2OI0I ; iFreeland fc Wilcots! 



l201lT ). higher than the density found from absorption-line 
measurements. These measurements rely on a set of assumed 
physical parameters including viewing angle, AGN proper 
motion, and kinematic luminosity. These assumptions are 
largely independent of those required for UV and X-ray 
density estimates, making this analysis complementary to 
existing methods. Bent-double radio sources also probe the 
density of the entire IGM, rather than just one temperature 
phase. 

Ram pressure resulting from the proper motion of these 
double-lobed radio sources through the IGM sweeps back 
the bipolar jets, producing the distinctiv e bent-double radi o 
tails no ted in many r adio observations (|Milev et a.1.1 1 1972T ) . 
Prior to lBurns et al.l (|l987f ) . it was believed that the condi- 
tions on ambient IGM density and galaxy velocities neces- 
sary to produce these bent-double radio tails could be found 
only in large, rich clusters of galaxies. However, surveys 
have found a significant number of bent-double sources in 
lower mass galaxy groups dVenkatesan et al .1 1994 1 Doe et al.l 
ll995l ; lBlanton et aUl200lf ). lEkerd (|l978h was among the first 
to identify bent-double radio sources as a possible density 
probe of the IGM, observing IGM densities in the range 
n ~ 3x lCT 4 -6x 10~ 3 cm" 3 in the galaxy groups NGC 6109 
and NGC 6137. These early estimates have been corrob- 
orated by more recent observations of head-tail sources 
in galaxy groups and po or clusters (jFreeland et al.l 120081 ; 
IFreeland fc Wilcotsl 1201 ll ) . If these numbers are represen- 
tative of all galaxy groups, a significant fraction of the local 
universe's baryon content would reside within the IGM. 

Analytic modeling of bent-double radio sources has 
found that the radius of curvature at the most bent part 
of the jet can be described by an equation of the form 



R 
h 



It 

Pra, 



(1) 



which balances internal and external pressure gradients as 
a ratio of the radius of curvature R and t he scale height h 
acros s which the pressure d i fference acts | Bcgclma n et al. 
19791; iJones fc Owenl 1 19791: iBurns fc Owen! Il98fj; lO'Deal 



19851 ). In iBegelman et all l| 19791 '), h is taken to be the di- 
ameter of the jet, whereas in Jones & Owen, h is the scale 
height of the ISM within the host galaxy. We use the diam- 
eter of the jet for h. We use this relationship as an analytic 
model to compare to our simulation results. 

We carry out a series of numerical simulations of bent- 
double radio sources to determine values for R and h in 
terms of initial model parameters. In §2, we describe the 



numerical methods and initial conditions used in our sim- 
ulations. In §3, we describe the results of our simulations 
and our methods for measuring R and h. In §4, we use our 
simulations to quantify the errors on density estimates from 
observations of bent-double radio sources. We also develop a 
formula to estimate the kinetic luminosity of observed jets, 
and we make predictions for the X-ray detectability of radio 
sources. In §5, we summarize our results. 



2 TECHNICAL DESCRIPTION 

2.1 Code description 

Simulations are carried out using the FLASH 2.4 hydro- 
dynamics code (|Frvxell et al.l 12000] ) . which is a modular, 
block-structure adaptive mesh code. It solves the Riemann 
problem on a three-dimensional Cartesian grid using the 
piecewise-parabolic method. Gas is modeled as having a uni- 
form adiabatic index of 7 = 5/3. 

2.2 The Jet Nozzle 

In order to simulate the injection of collimated, supersonic 
jets into the grid, we empl oy a numerical "no zzle", as first 
developed and described in I Heinz et al l (|2006l l: an internal 
inflow boundary of cylindrical shape placed at the location 
of the AGN, injecting fluid with a prescribed energy, mass, 
and momentum flux to match the parameters we choose for 
the jet. For reasons of numerical stability, we impose a slow 
lateral outflow with low mass flux in order to avoid complete 
evacuation of zones immediately adjacent to the nozzle due 
to the large velocity divergence at the nozzle. The injection 
of energy and mass due to this correction is negligible. 

We model unresolved dynamical instabilities near the 
base of the jet by imposing a random-walk jitter on jet axis 
confined to a 5° half-opening an gle. This is ne cessary to 
model the 'dentist's drill' effect of lScheuerl (|l982T l. 

We chose to inject the jet at an internal Mach number 
of 10 in most simulations. For computational feasibility, we 
chose a jet velocity of Vj e t = 3 x 10 9 cm/s for most simula- 
tions. The jet is turned on initially and continues to inject 
material for the entire length of the simulation. 



2.3 Initial Conditions 

O ur setup is s imilar to that used for X-ray binary jets 
in lYoon et alj (|201ll ). We place the jet nozzle in a mov- 
ing medium inside a simulation box large enough that the 
boundaries never affect the simulation (2.8 Mpc on a side). 
We keep the location of the AGN fixed in space, letting 
the IGM stream by at velocity — n ga i perpendicular to the 
jet axis. We can vary the velocity and density of the IGM, 
as well as the luminosity, velocity and internal Mach num- 
ber of the jet. In all cases the IGM pressure is set to 
2.76 x 10 -13 erg/cm 3 . A list of parameters for all simula- 
tions is given in table [T] 

The simulations were carried out on a staggered mesh 
grid as described in lYoon et al.l (|201ll ) in order to capture 
the large dynamic range required, and ensure that the nozzle 
diameter is resolved by 12 grid cells. For our normal scaling, 
the nozzle diameter is 2 kpc, with a maximum resolution 
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for the standard model of about 0.175 kpc near the jet noz- 
zle. The nozzle diameter for each simulation, dj, is listed 
in table [T] In all cases, the maximum resolution is set such 
that the nozzle is resolved by 12 grid cells. Short duration 
test simulations at higher resolution were carried out and 
produced similar results. 

One simulations was for a conical jet with an initial 
opening angle, rather than a purely parallel inflow. This 
model, .015E_theta20, has identical parameters to model 
.015E, except that the nozzle is 4 times smaller (and the 
maximum resolution 4 times greater) and the inflowing ma- 
terial is spread over a 20° half-opening angle. This model 
was run as a direct comparison to model .015E, based on 
the analytic model in section [4~T1 We use cylindrical (paral- 
lel) jet injection for our other models because this allows us 
to use a much lower resolution. 



3 SIMULATION RESULTS 

When the jet initially turns on, it drives a shock into the 
IGM, creating a rapidly expanding cocoon. The expansion 
velocity of the cocoon is initially faster than w ga i, so the AGN 
remains inside the cocoon and the jet stays fairly straight. As 
the cocoon expands it decelerates, but u ga i remains constant. 
Eventually, the AGN moves outside the initial cocoon and 
a bow shock is created around the AGN jets. The pressure 
gradient from this shock bends the jets backwards, creat- 
ing the characteristic curved structure. At some point along 
the jet, the flow becomes unstable and breaks up, creating 
bright radio lobes at the ends of the jets. The swept-back jet 
material creates a long tail of diffuse jet material reaching 
back to the cocoon created around the initial position of the 
AGN. 

We create synthetic radio images by assuming that 
the jet material, followed with a tracer fluid, consists of a 
hot plasma of relativistic particles with magnetic fields in 
equipartition with thermal pressure. Energy lose from syn- 
chrotron cooling is not included. 

The left panel of Fig. [1] shows a synthetic radio image 
for simulation .25E at 50 Myr. In this image the AGN is 
moving relative to the IGM from right to left. The narrow 
curved jets and bright radio lobes can be clearly seen. The 
tail and cocoon produce diffuse radio emission with a surface 
brightness 1 to 2 orders of magnitude fainter than the jet and 
lobes. This extended emission is typically not seen in bent- 
double radio sources in jets. It is possible that this emission 
is resolved out in existing observations or that synchrotron 
cooling, which is not included in making our images, makes 
the emission at 1.4 GHz too faint to observe. Observations 
at lower resolution and/or lower frequencies may reveal the 
extended tail. It is pos sible that at least p a rt of the tail 
is seen in source S7 in iFreeland fc Wilcotd (|201 J ), which 
appears to have extremely thick jets relative to the radius 
of curvature. However, this is not the only explanation for 
this source's appearance (see section PO)) ■ 

The length of the tail is v ga i x £agn , where *agn is 
the amount of time the AGN has been active. If there is 
an estimate for the AGN host galaxy velocity, finding the 
length of the tail would allow the amount of time the AGN 
has been active to be determined. Note that in Fig. [1] the 
length of the tail (left to right) is less than the width of 



the tail (top to bottom), because the jets initially expand 
outwards faster than v ga i. If the AGN remains active, the 
tail will continue to lengthen and eventually become longer 
than it is wide. 



3.1 Radius of Curvature 

The right panel of Fig. [T]shows the same image with a circle 
overlaid with R — 37 kpc, the best fit radius of curvature 
for this image. The circle traces both the upper and lower 
jet and passes through the bright radio lobes where each jet 
breaks up. Note that the upper lobe is significantly brighter 
than the lower lobe at this time. The brightness of the lobes 
and the curvature of the two jets varies with time due to 
instabilities in the jet propagation and the small random 
changes in jet direction that we introduce. 

The radius of curvature is determined by an automated 
fitting procedure. We start with a synthetic radio image and 
consider only the region 28 kpc above and below of the AGN, 
the approximate extent of the jets. Because the upper and 
lower jets and lobes can differ significantly in brightness, we 
find the maximum brightness of any point in each half of the 
image and exclude points with less than 10% of this bright- 
ness. For each slice along the jet, we then determine the 
horizontal location of the jet by taking an intensity weighted 
average of the radio emission in that slice. We then fit a cir- 
cle to the jet locations weighted by the total radio intensity 
of the points in each slice. 

To characterize how much fluctuations in the jet affect 
the radius of curvature, we use the above procedure to mea- 
sure the radius of curvature at a series of times (after the 
jet curvature is well established) and compare the results. 
Fig. [2] plots measured values of R for simulation .25E from 
40 Myr to 220 Myr. The error bars represent the 1-cr error 
on the value of R at each time, typically about 10%. The 
values of R are very consistent, with a scatter of about 25%. 
The best-fit value of R for simulation .25E is 35.5 ±7.9 kpc, 
where the error takes into account both the error in indi- 
vidual fits and variations between different snapshots. The 
same procedure is applied to all of our simulations, with the 
results listed as R in table [TJ 



3.2 Jet Thickness 

A similar process is used to find the average jet thickness, 
listed as h in table Q] The thickness typically has a very 
small variation, ~ 5%, and the accuracy is limited by the 
resolution of our simulations. Because we are interested in 
the "real" value of the jet thickness, we determine it using 
the raw simulation data rather than the synthetic radio data. 
We take a region of ±24 kpc from the AGN and define points 
as being in the jet if at least 80% of the material at that point 
was injected by the nozzle and it is moving with at least 80% 
of the initial jet velocity. This excludes all material in the 
lobes and the cocoon around the jet. We then determine the 
jet thickness perpendicular to both the jet direction and the 
motion of the AGN, henceforth the 'z' direction, as this is 
the front along which the ram pressure acts. For each slice 
through the jet we take the average thickness for points in 
the jet in the 'z' direction. We then take the average of this 
thickness in all slices to get a value for h at that time, with 
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Table 1. Model Data 
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.07 


2.0 


lE_4n 


1.0 


4 x 10" 


3 


1000 


0.1 


10 


35.5 ± 12.2 


2.09 ± 
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a Nozzle diameter at injection 

b This simulations has a half-opening angle of 20° at injection, rather than a parallel inflow 



the standard deviation of the average being the uncertainty 
at that time. To get the values of h in table [1] we then find 
the best-fit value of the thickness for all times considered, 
and the error takes into account both the error at individual 
times and the variation with time. 

The overall result is that the radius measured at a par- 
ticular time is accurate to within an error of 25% of the 
'true' value for that particular combination of parameters. 
The thickness is more consistent, with typically only a 5% 
variation. The variation in R sets a lower limit on the ac- 
curacy of the IGM density estimated from observations of 
bent-double radio sources. 



4 DISCUSSION 

4.1 Analytic Fit for Radius of Curvature 

For an AGN moving relative to the IGM, the jets will curve 
due to the ram pressure of the incoming IGM. At the point 
of maximum curvature, the ratio of the radius of curvature 
to the thickness of the jet will equal the ratio of the ram pres- 
sure of the jet to the external ram pressure dBegelman et all 
1 19791 ; I Jones fc Owenll 19791 ; iBurns fc Owenlll98fj ). i.e.: 



jet 



R _ 

h piGMVga] 

The ram pressure of the jet will be 



-Piet — 



-jet 



fft 2 «jet 



(2) 



(3) 



Rearranging eqn. [2] and eqn. Ogives a formula for R in 
terms of h: 



R 



Ljet 



f hVjetpiGMV^ 



(4) 



However, the thickness of the jet is not constant and will 
be determined by the external pressure. Initially the lateral 
expansion of the jet will be ballistic because the component 
of the ram pressure perpendicular to the jet will be higher 
than the external pressure. However, as the thickness in- 
creases the internal pressure will drop until it is the same 



order as the cocoon (i.e. external) pressure. At this point, 
a recollimation shock will be driven into the jet, setting the 
thickness and providing an internal pressure that is equal to 
that of the cocoon. Therefore, the thickness will be set such 
that the external pressure balances the perpendicular ram 
pressure. For a jet with an initial half-opening angle off), the 
average ram pressure of the jet perpendicular to direction of 
motion will be Pj e t,± — ^-Pj c tsin 2 #. Balancing this against 
the external pressure gives 



piGMWgai = --Pjct sin 2 



1 ijet . 2 a 

sin 



2 -j- 2f/l2«j et 



Solving for h, we find 



h = 



1/2 



V2 V Z«j<=tpIGMfg al 

Using this value for h in eqn. [4] gives us 



(5) 



(6) 



7? = 



2h 
sin 2 6> 



V2_ 
sin 9 



r UjctPlGMW ga l 



1/2 



(7) 



For our simulations with a cylindrical jet, we fix the 
nozzle size to h set by 



4 kpc x 



rc-aWgai.iooo^t.o.i 



1/2 



(8) 



where L 44 = L jct /10 44 erg/s, rt_ 3 
fgai.iooo = Ugai/1000 km/s, and Wj ot ,o.i 



" = n 1G M/W~ cm- 3 , 
= Vj et /(0.1c). This 
is the equivalent width of a jet with an initial opening an- 
gle of 6 = 20°. From eqn. this would predict a value of 
h/R = 1/17. 

We run one model, .015E_theta20, with a small nozzle 
size and a 20° initial opening angle. The measured values 
of R and h in talbe[T]are 8.6 ± 0.8 kpc and 0.44 ± 0.2 kpc, 
respectively, very close to the prdicted values of 8.5 kpc and 
0.49 kpc given by equations [7] and [U The results are also 
very similar to model .015E, which has the same setup but 
with a cylindrical jet input and a 0.5 kpc nozzle size. 

Using the measured values of R and h in table [T] for all 
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of our parallel jet models, we find that the value of h stays 
very close to the nozzle size and the ratio of h/R is about 
1/18, close the our analytic estimate of 1/17. 

Using the measured values of R in table [1] we find that 
the value of R in terms of our model parameters is 



R = 5 



72 kpc x 



1/2 



n -3Vl,1000 V jet,0.1 



1/2 



(9) 



Figure [3] plots the radius of curvature vs. jet luminosity 
for models .015E, .062E, .25E and IE. These models differ 
only in jet luminosity, with all other parameters the same. 
Error bars represent the overall error in measuring R, as 
discussed in section |3~T1 The solid line is a plot of eqn. [5] for 
radius vs. luminosity and passes well within the error bars 
of all the four data points. 

Equation [7] can also be used to find the kinetic lu- 
minosity of observed bent-double radio sources based on 
their measured rad i us of curvature and jet pressure. In 
iFreeland k, Wilcotsl (|201ll ). Pj e t is assumed to be the mini- 
mum synchrotron pressure 



Pjet — Pmin — 



(2^) 3/7 (J^j [c^L^il + k^V)- 1 ] 4 ' 7 (10) 

where C12 is a constant t hat depends on t he spectral in- 
dex and frequency cutoffs ljPacholczvklll970l ). k is the ratio 
of relativistic proton to relativistic electron energy, (f> is the 
volume filling factor, V is the source volume, and L ra( j is the 
radio luminosity of the jet at the point where the pressure 
is measured. Pj c t is measured in several slices along the jet, 
with the volume assumed to be proportional to the square 
of measured jet thickness, which is limited by the beam size 
of the observations. Therefore, P m in oc h~ 8 ^ 7 but is inde- 
pendent of the radius of curvature R and the length of the 
jet. 

Rearranging eqn.0 we find a jet luminosity of 



Ljet = -^R 2 sin 2 OpiGMV^iV 



jet 



= 3.73-Rkpc mGM,-3Ugai,iooo/3jet x 10 42 erg>i) 

where Pk P c is R in kpc and /3j c t = fjet/c. In terms of pressure 
Ljct is 



2.24/ig pc Pjet,-ii/3 jct x 10 42 erg 



(12) 



where h^ pc is h in kpc. 

The jet kinetic luminosity calculated using this formula 
for the sources and measur ed values of R, h a n d the syn- 
chrotron pressure P m m in IFreeland fc Wilcotsl (|201ll ) are 
listed in table [2] The values in table [2] are upper limits on 
the luminosity made with the assumptions that the observed 
value of h is the true jet thickness and that the jet velocity is 
Vj ct = c. If the jet is narrower or slower, the luminosities will 
be smaller. Note that the formula for Lj c t is independent of 



the AGN velocity, so values can be found even for sources 
with unconstrained velocities. The sources are all of order 
10 45 erg/s and the variation between sources is smaller than 
the variation in the total radio power at 1440 MHz (L1440). 
Pmin is proportional to L^fJ d , which is the total radio emis- 
sion at the point where P m i n is measured. I/ ra d therefore de- 
pends on the 1440 MHz emission at that point and the model 
of the synchrotron spectrum used, but is not directly depen- 
dent on the total 1440 MHz emission of the entire source. 
Pmin scales with jet thickness as h so -Ljct oc 

Note that we assume Pj e t comes only from synchrotron 
pressure. If the true jet pressure is higher, the estimates of 
both niGM and Lj ot in table[2]will be higher. Adding invisible 
components to the momentum flux of the jet only acts to 
increase the required IGM density. On a similar note, in our 
simulations very little external material becomes entrained 
in the jets. However, even if a significant amount of mass is 
entrained it will not change the momentum flux of the jet. 
Therefore, the radius of curvature and inferred IGM density 
should not affected. 



4.2 Effects of Observational Resolution 

Although the jets are well resolved in our simulations, 
they are typically unresolved in radio observations. From 
eqn. [2] and eqn. 1101 the density derived from observations 
will scale with jet thickness and radius of curvature as 
n-iGM oc (/i/P) _1//7 P _8//7 . The values of h/R in our simula- 
tions (table[TJ range from about 1/13 to 1/21 with a typical 
value of about 1/18. The observed ratio (table[2} ranges from 
h/R = 1/14 for S3 to h/R = 1/0.83 for S7. Although the 
density scales weakly with jet thickness, if we assume a real 
value of h/R — 1/17, corresponding to an initial jet opening 
angle of 6 — 20°, the densities for observed sources would 
be between 3% (S3) and 54% (S7) higher. If a = 5° inital 
opening angle is assumed, the correction would be between 
52% (S3) and 128% (S7) For an under-resolved jet, the ob- 
served thickness will always be too high, and therefore the 
density derived will always be lower than the actual value. 
This is generally a fairly small correction, about 25% for a 
typical source, but can be up to 50% or more for sources 
with ftobs ~ R- 

Inadequate resolution can also affect the measurement 
of the radius of curvature. To characterize this, we apply a 
Gaussian smoothing filter to our radio images of simulation 
.25E to simulate radio beam FWHM sizes from 1 kpc to 
32 kpc, and then used our fitting routine to find the radius 
of curvature. The results are plotted in Fig. [3] For this sim- 
ulation, the actual radius of curvature is R = 36 kpc and 
the jet thickness is h = 2 kpc. For a marginally resolved 
jet (beam size ^ 2 kpc) the measured value of R does not 
change. For larger beam sizes, R is over-estimated by about 
10% to 25%. As the beam size becomes comparable to the 
radius of curvature, the fit for R becomes very poor (error 
comparable to R) . Even in this case, however, the measured 
value of R is systematically larger than the real value. An 
over-estimate of R will lead to an underestimate of the den- 
sity derived from eqn. [2] typically about 20% for sources 
with unresolved jets. 

Combining the effects of over-estimating the jet thick- 
ness and over-estimating the radius of curvature, we find 
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Table 2. Jet Kinetic Luminosity 



Source ID 


a.b 

V ' 
gal 

(km /s J 


pa 

10 ^ ^rg cm ^ 


IGM 
f cm ^ "1 


h a 


f>a 


T a 

^1440 

(W Hz -1 ) 


( erp" /si 


SI 
S2 
S3 
S4 
S5 
S6 
S7 


430±J™ 
570 ± 60 d 
745+109 
n 50 +210 

Unconstrained 
Unconstrained 
850l 120 


0.9 ± 0.2 
0.6 ± 0.2 
1.4 ±0.6 
1.7 ± 0.3 
0.4 ±0.1 
0.6 ±0.1 
1.4 ±0.3 


3± 2 x 10" 3 
5±4x 10~ 4 
2± 1 x 10~ 4 
5± 2 x 10" 4 
(23 ± 7)/v 2 
(52 ± 16)/v 2 
2 ± 1 x lO" 3 


23 ±1 
30 ± 4.5 
10 ±0.6 
22 ±0.7 
38 ±3.8 
19 ±3.1 
22 ± 3.6 


42 ± 6 
104 ±9 
141 ± 19 
89 ± 7 
220 ± 11 
69 ±4 
18 ±4 


1.6 x 10 25 
1.88 x 10 25 
3 x 10 23 
9.5 x 10 24 
6.4 x 10 24 
9 x 10 23 
2 x 10 24 


1.1 ± 0.2 x 10 45 

1.2 ±0.5 x 10 45 
3.1 ± 1.4 x 10 44 
1.5 ±0.3 x 10 45 

1.3 ± 0.4 x 10 45 
4.8 ± 1.4 x 10 44 
1.5 ±0.5 x 10 45 



a from lFreeland fc Wilcotsl |201ll) . 

b v gal here is \/3 times the group velocity dispersion, which is listed as v ga \ in Table 1 of Frecland fc Wilcotsl (l201lh . 
c Upper limit, assuming tij et = c 

d For S2, velocity is based on the difference in redshift between the source galaxy and the group, not the velocity dispersion. 



that for a typical source in iFreeland fe Wilcotsl I^OllT l (re- 
solved source, unresolved jet, h/r ~ 1/4) the density cal- 
culated is low by about 50%, assuming a true value of 
h/R of 1/17. This ranges from no correction for source S3 
(marginally resolved jet) to about 85% low for source S7 
(beam size ~ R). 

For estimates of the jet luminosity, the error due to 
resolution can be significantly larger. From eqn. 1121 we find 
that I/j c t oc h Pjet oc h 6/r . There is no dependance on R, 
but luminosity is very sensitive to h. For example, assuming 
a real value of h/R of 1/17, the derived luminosity in table[2] 
would be reduced by a factor ranging from 1.2 (S3) to 13 
(S7). 

4.3 Effect of Viewing Angles 

So far, we have produced synthetic images of our simulations 
assuming that both the direction of jet propagation and the 
direction of motion of the AGN relative to the IGM are 
perpendicular to the observer's line of sight. In reality, for 
observed bent-double radio sources there will be unknown 
angles between the jet direction and the observer and the 
direction of motion and the observer, both of which will 
affect the measurement of the radius of curvature. The an- 
gle between the jet and the direction of motion will change 
where the point of maximum curvature is along the jet, but 
will not change the r adius of curvature at that point (e.g. 
iBegelman et alJll979l ). 

To quantify how much error viewing angle introduces, 
we take the data from simulation .25E, rotate it by various 
angles, then use it to create synthetic radio images. We then 
measure the radius of curvature in each of these images by 
the procedure in sect. 13. ll Fig. [5] plots the measured radius 
of curvature for an angle between the jet direction and the 
observer between 0° and 90°. In all cases the direction of 
motion of the AGN is perpendicular to the observer's line of 
sight. The solid line in Fig. [5]is R = R(g=o) x cos(t9), which 
is the expected apparent radius of curvature of a simple 
rotated circle. The measured value of R follows this line 
fairly well out to about 60°, where the apparent radius is 
half the real radius. Beyond this, the measured radius does 
not get any smaller. This is because the radio lobes at the 
ends of the jets begin to overlap from the observer's point 
of view, and their size dominates the fitting routine. This is 



one possible explanation f or the appearance of source S7 in 
IFreeland fc Wilcotsl (|201lh . which appears to have a small 
radius of curvature but with very thick jets. At 90°, the 
fitting routine fails. However, it is unlikely that a source 
with the jet aimed almost directly at the observer would be 
classified as a bent-double radio source. 

Fig. [6] plots the measured radius of curvature for an 
angle between the direction motion of the AGN and the 
observer between 0° and 90°. In all cases the jet direction 
is perpendicular to the observer. The solid line is Fig. [6] is 
R — 7?(0 =O )/ cos(^), which is the expected apparent radius 
of curvature of a simple circle rotated in the same manner. 
The measured value of R follows this line out to about 60°, 
where the apparent radius is double the real radius. Beyond 
this, the source is so inclined that the jet appears the be 
approximately straight and the fitting routine fails. 

In practice, sources will be rotated around both axes, 
and orientation effects tend to cancel each other out some- 
what. Assuming sources rotated beyond 60° in either direc- 
tion will not be classified as bent-double radio sources, the 
unknown viewing angle still introduces a large uncertainty 
in the measurement of R, and therefore in the values of uigm 
and I/jct- On average, the error in the measurement of R due 
to viewing angle should be between 50% low to 30% high 
(1-c), but for an individual source could be up to a factor 
of 2 (100% under-estimated to 50% over-estimated). 

4.4 X-ray Detectability 

In general, the gas in galaxy groups is too cool and too low 
density to be seen in X-ray observations. However, as an 
AGN moves through the IGM, a cocoon of shocked material 
develops around the jets. The shock heats and compresses 
the IGM, boosting the X-ray emission. The parameters that 
affect the X-ray brightness are the velocity of the AGN, 
Vgal, which determines the temperature of the shocked gas, 
and the density of the IGM, tiigm, which will determine the 
emissivity at a given temperature. 

To determine if bent-double radio sources in groups 
would be detectable in X-rays, we used the XIM tool 
ijHeinz fc Briiggenl [2009) to model 100 ks Chandra obser- 
vations of four of our simulations, models .015E_.25vel, 
.062E_.5vel, .25E, and lE_4n, placed at a redshift of z = 0.1. 
Fig. [7] shows synthetic Chandra images. In all figures we as- 



© 0000 RAS, MNRAS 000, 000-000 



Simulations of bent-double radio sources in galaxy groups 7 



sume a background IGM temperature of 5 x 10 5 K and an 
IGM metallicity of Z — 0.3 solar. All four simulations have 
about the same radius of curvature of R « 36 kpc. The first 3 
simulations have the same IGM density with AGN velocities 
of Ugai = 250 km/s (upper left), 500 km/s (upper right) and 
1000 km/s (lower left). The last two simulations have the 
same AGN velocity (u ga i = 1000 km/s) but different densi- 
ties of niGM = 10~ 3 cm~ 3 (lower left) and 4 x 10~ 3 cm' 3 
(lower right). For these images, we assume there is no point- 
source emission from the AGN. Even if there is point-source 
emission, however, the high spatial resolution provided by 
Chandra would be able to remove this contribution. 

Except in the lowest velocity case, there is significant 
X-ray emission from the region of the AGN jet and extended 
tail. We would expect emission to be strongest at the leading 
edge of the jet, where we are seeing the strongest part of the 
shock edge-on. However, because the edge-on shock is very 
thin and significant smoothing is needed to bring out the 
X-ray emission (10 arcsec in these images), the leading edge 
does not stand out. Instead, X-ray emission is spread across 
the region affected by the AGN and dominated by face-on 
shocks. 

In Fig. [8l we plot the mean surface brightness between 
—50 and +50 arcsec of the AGN along the y-axis (parallel 
to the jet) in Fig. [7) In the worst-case presented here (model 
.015E_.25vel, upper left panel), the surface brightness is 
about 0.02 counts/arcsec 2 for a 100 ks exposure, about 
twice the background level. As the AGN velocity increases, 
the surface brightness increases as roughly vZ^ , reach- 
ing about 0.045 counts/arcsec 2 for model .25E. For model 
lE_4n (lower right panel), which has the same velocity as 
model .25E, the brightness increases to 0.3 counts/arcsec 2 , 
a scaling of about ^igm • 

If these models were placed at a higher redshift, the an- 
gular size of the X-ray source would decrease, but the surface 
brightness would remain the same. For a large radius of cur- 
vature, the thickness of the shocked material, and therefore 
the surface brightness, will increase proportional to R. 

X-ray observations of bent-double radio sources in 
groups would place important constraints on the IGM prop- 
erties. The X-ray surface brightness, along with an estimate 
of Dgai from velocity dispersion in the group, would allow 
tiigm to be calculated independent of the radio observations. 
Density calculated this way would also be less sensitive to 
the AGN velocity, scaling as about wigm oc v^ ' 5 rather 
than niGM oc as in eqn. [2] The different scaling also 
means that in cases where w ga i is unconstrained, X-ray and 
radio data can be combined to find both the IGM density 
and AGN velocity. This can also be used to refine the vales 
of w ga i and wigm and constrain the viewing angle in cases 
where there is an estimate for the AGN velocity. Because the 
X-ray emission follows the tail of extended radio emission, 
measuring the length of the X-ray emitting region would 
provide an age estimate for the source. 

The general formula for the X-ray surface brightness is 

S~0.1x (Z + 0.1) x Wgjfiooo x ri-l x (-R/36 kpc) 

counts/arcsec 2 (13) 

for a 100 ks observation, where Z is the metallicity of the 
IGM relative to solar. For sources SI and S2 in table[2]there 
are Chandra observations (|Freeland et alj|2008h which found 



a total of 80 ± 35 and 94 ± 26 counts above the background 
in a 35.17 ks and 47.19 ks exposure, respectively. Assuming 
the total area on the sky of these sources is about 4i? 2 (the 
actual size of the emitting region is unknown), the total 
counts predicted would be very roughly 30 counts each for 
Z — 0.3. This is lower than the observed counts, but current 
observations are not good enough to distinguish counts from 
the region of the bent-double from rest of the IGM, so they 
are still consistent. 

For particularly bright sources, there may be enough 
photons to determine the temperature of the shock. The 
temperature will scale roughly as T oc w 2 al , so if the tem- 
perature can be found it would allow w ga i and tiigm to 
be estimated from the X-ray data alone. With additional 
constraints from radio and velocity dispersion observations, 
it would then be possible to find the compression ratio of 
the shock, and from that the temperature of the unshocked 
intra- group medium. 



5 CONCLUSIONS 

Our simulations are able to closely match the appearance 
of observed bent-double radio sources with narrow, curved 
jets ending in bright radio lobes. We also predict that there 
should be an extended tail of radio emission that may be 
observable at low resolution and/or low frequencies. The 
length of this tail would allow the age of the AGN to be 
determined. 

From our simulations, we derive a formula for the radius 
of curvature (eqn. [5} in terms of IGM density, AGN velocity, 
jet luminosity and jet velocity. From this we are able to 
calculate the kinetic jet luminosity (eqn. I12[l for observed 
bent-double radio sources, and find that Lj et is typically 
around 10 45 erg/s, assuming that Vj e t — c and that the jets 
are really as thick as their observed values (see table [2}. 
The luminosities should be considered upper limits, as they 
will be lower if either the jets are slower or the jets are 
intrinsically thinner. This formula is independent of v sa i and 
therefore can be used to find Lj et even when the velocity of 
the AGN is unknown. 

A lack of resolution in radio observations leads to a sys- 
tematic under-estimate of the IGM density for two reasons: 
the jet is unresolved and the radius of curvature is over- 
estimated. In our simulations, use initial conditions that 
produce ratio for jet thickness to radius of curvature of 
h/R ~ 1/17, equivalent to a 20° initial opening angle of the 
jet. For observed sources (limit ed by resolution), th i s rati o 
ranged from 1/14 to 1/0.83 in IFreeland fc Wilcotsl (|201ll ). 
Density scales as (/i/i?) _1//7 , so this leads to an under- 
estimate of the IGM density of up to 50%. Inadequate reso- 
lution also leads to an over-estimate of the radius of curva- 
ture (and corresponding under-estimate of density), proba- 
bly by abo ut 20%. Overall, IGM dens ity estimates for typical 
sources m IFreeland fc Wilcotsl i|201ll ) are low by about 50% 
due to resolution effects. This ranges from no correction for 
source S3 (assuming a marginally resolved jet) to about 85% 
low for source S7 (beam size ~ R). 

The largest source of error, however, comes from the 
unknown angles between the observer, the jet direction and 
direction of motion of the AGN. This can lead to either an 
over- or under-estimate of the true radius of curvature. This 
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leads to an uncertainty of up to a factor of 2 in the esti- 
mate of the IGM density, with a typical (1-cr) error or about 
50%. This is comparable to the error from all observational 
uncertainties (dominated by the uncertainty in the AGN ve- 
locity). Even this uncertainty, however, does not change the 
conclusion that bending can be used to diagnose pigm in a 
statistical sense. 

Finally, we have modeled the X-ray emission of bent- 
double radio sources and predict that they should be de- 
tectable in Chandra X-ray observations. Although the IGM 
in groups is generally too cool and diffuse to be seen in X- 
ray observations, the shocks around the jet and tail in our 
simulations compress and heat the IGM, potentially making 
the entire region affected by the AGN jets detectable in X- 
rays. The X-ray surface brightness scales as approximately 
S oc wJo^DgJ 5 . Sources in fairly dense environments and 
with fairly large angular size should be detectable in moder- 
ate (~ 100 ks) observations. Count rates fro m existing short 
X-ray observations of sources SI and S2 in iFreeland et all 
(2008) are consistent with our predictions. Further X-ray ob- 
servations would provide complimentary constraints on the 
IGM density and AGN velocity to radio observations. 

For particularly bright sources, it may be possible to 
obtain a measure of the temperature of the shock from X-ray 
observations. This would provide an independent measure of 
the AGN velocity relative to the IGM. 

Future X-ray and radio observations will be able to 
place better constraints on the IGM density, AGN velocity 
and AGN age. With very complete observations it should 
also be possible to constrain the temperature of the IGM 
and the orientation of the jets and direction of AGN mo- 
tion. 
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Figure 1. Left: Simulated radio emission from model .25E at 50 Myr (log scale). Color bar shows intensity in mjy/arcsec 2 and brightness 
temperature at 1440 MHz. The image is 170 kpc on a side. In addition to the two bright radio jets, there is diffuse synchrotron emission 
filling the cocoon created by the jets. The surface brightness of this emission is ss 10 — 100 times fainter than the jets. Right: Same as 
left panel, but overlaid with the best fit radius of curvature (white circle). The radius is 37 kpc. 
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Figure 2. Best fit radius of curvature at different times for model .25E. The radius varies about 15% between measurements. The 
overall mean radius is 35.5 kpc with an error of 7.9 kpc. 



© 0000 RAS, MNRAS 000, 000-000 



10 B. J. Morsony et al. 



1 00 



o 

Q_ 



80 



<>/ 



CD 



D 
> 



o 



TD 
D 



60 



40 



20 



,42 



■0 



4.5 



10 fl 



L (erg/s) 



Figure 3. Radius of curvature vs. jet luminosity (diamonds) for models .015E, .062E, .25E and IE. These models differ only in jet 
luminosity, with all other parameters the same. Error bars represent the overall error in measuring R. Solid line is a plot of eqn. [9] for 
radius vs. luminosity. 
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Figure 4. Measured radius of curvature vs. beam size for model .25E. Error bars represent the overall error in measuring R. As the 
beam size increases, the measured R initial remains fairly constant, but then begins to increase, along with the error, as the beam size 
approaches R. Once the beam size is approximately equal to R, the radius cannot be reliably measured. 



© 0000 RAS, MNRAS 000, 000-000 



Simulations of bent-double radio sources in galaxy groups 11 



60 



o 




100 



Angle of jet to observe (degrees) 



Figure 5. Measured radius of curvature vs. angle of jet direction relative to the observer (0° is perpendicular) for simulation .25E. Solid 
line is R = R(g- ^ X cos(0), the expected geometric correction. The measured R follows this line fairly well out to about 60°, beyond 
which the radius of curvature is poorly fit. 
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Figure 6. Measured radius of curvature vs. angle of direction of motion relative to the observer (0° is perpendicular) for simulation 
.25E. Solid line is R = R^—^/ cos(</>), the expected geometric correction. The measured R follows this line fairly well out to about 60°, 
beyond which the radius of curvature is poorly fit. 
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Figure 7. Simulated Chandra X-ray observations for a lOOks of models .015E_.25vel (upper left), .062E_.5vel (upper right), .25E (bottom 
left) and lE_4n (lower right) at a redshift of z = 0.1. All images are the total counts from 0.3 to 3 keV and are smoothed over 10 arcsec. 
The radius of curvature of the jet in each model is about 36 kpc (20 arcsec). Color bars show the scale of each image in counts/arcsec 2 , 
scaled to 70% of the maximum brightness in each image. Orange contours mark the outer extent of radio emission. 
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Figure 8. Plots of the mean surface brightness between y = —50 and y = +50 arcsec in each images in Fig.[7]vs. x position. Even in the 
worst case (model .015E_.25vel, upper left), the mean surface brightness is about .02 counts/arcsec 2 , about twice the background level. 
This increases roughly as v®£p to about .045 counts/arcsec 2 for model ,25E. For model lE_4n, which has the same velocity as model 
.25E, the brightness increases to 0.3 counts/arcsec 2 , scaling as about ^jq m - 
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